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SUMMARY 

This work investigates the effects of a force induced by the labyrinth 
seal on the stability of rotor systems and the factors of the seal which affect 
the stability. In the analysis, it is assumed that the fluid in the seal is 
steady and that the rotor is set vertically in order to avoid the effects of 
gravity force. The force induced by the seal is expressed in terms proportion- 
the velocity and displacement of the rotor and is deduced to that expres- 
sion for the oil-film force in journal bearings. That force is taken into ac- 
count in the equations of motion; then the stability of the system is discussed 
by energy concept. 

The force induced by the labyrinth seal always makes the rotor system un- 
stable, and the tendency is remarkable when seal leakages are small. The reso- 
nance point of the rotor system is also affected by the labyrinth seal; that is, 
the resonance point of the rotor system is removed by the seal leakages. The 
flow pattern in the labyrinth seal was investigated experimentally, and the 
force induced by the labyrinth seal was measured by using a water— tunnel experi- 
mental system which was designed to measure the labyrinth seal force by using 
the similarity between gas and liquid flow theory. 

INTRODUCTION 

After the oil shock, high-performance turbines and compressors are required 
in order to save energy. For this purpose, designers would like to minimize leak- 
age from labyrinth seals, so they design the clearances of the labyrinth seal 
to be small. However, if the clearances are small, self-excited rotor vibra- 
tions are caused by the flow forces of the working fluid. The origins of the 
exciting forces are at present only partially known as a steam whirl excitation. 
So it is not enough to evaluate these forces in order to design the labyrinth 
seals for compressors and turbines. Thus the analysis of labyrinth seals and 
the materials for design are strongly required by the designer of turbines and 
compressors. This paper is devoted to a basic analysis of the fluid force due 
to labyrinth seals. 

First, the fundamental equation proposed by Kostyuk (refs. 1,2) is extended 
in order to consider the effect of the variation of gland cross section. For 
the analysis, the fundamental equation is rewritten to ordinary differential 
equations by using the finite difference method. Then spring and damping co- 
efficients of the labyrinth seals are calculated for selected models from the 
fundamental equation and perturbation from the steady state. The flow rate and 
pressure, etc., in the steady state are also calculated. Then the stability of 
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the rotor system is discussed in terms of these coefficients by using the con- 
cept of energy. Furthermore, experiments were executed to observe the flow pat- 
tern in the gland and to study the characteristics of the flow-induced forces 
in the labyrinth seals. 

FUNDAMENTAL EQUATION FOR GAS FLOW IN LABYRINTH SEAL 

For the derivation of the equations, the following conditions are assumed; 

(1) Fluid in the labyrinth seal is assumed to be gas, and its behavior is 
assumed to be ideal. 

(2) Temperature of the fluid in the labyrinth seal is assumed to be con- 
stant . 

(3) Cross-sectional area of the seal gland is assumed to be constant in 
spite of the deflection of rotor, and only the time derivative of the cross- 
sectional area is considered. 

(4) Change of flow state in the gland is assumed to be isentropic. 

Thus from the illustrations of figure 1 the fundamental equations with respect 
to the flow rate and pressure are as follows; 
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Equations (1) to (7) are nonlinear partial differential equations; so for the 
analysis they must be rewritten to linear partial differential equations by 
using the perturbation terms from the steady state. Therefore, pressure, axial 
flow rate, and peripheral velocity in the steady state should be obtained. 

ANALYSIS OF FLOW IN STEADY STATE 

As the flow in the seal is steady, all state variables are constant; there- 
fore time derivatives and space derivatives of the state variables are zero. 

Thus the fundamental equations become as follows: 

( 8 ) 


?utiCi - JiCi-1 (9) 

From these equations, state variables in the equilibrium condition are obtained 
by using the iterative method as shown in figure 2. For this calculation, the 
pressure recovery factor rip is 
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and H is the angle between the rotating axis and the flow direction passing 
through the seal strip. 

The steady-state flow rate is given as 
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LINEARIZATION OF FUNDAMENTAL EQUATION 

For the linearization of the fundamental equation, the perturbations of 
pressure, peripheral velocity, and flow rate from those of the steady state are 
introduced as 


( 12 ) 
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where P*i, C Ai , and q Ai are pressure, peripheral velocity, and axial flow 
rate of steady state in the i*-^ 1 gland and hi> and Cl are the nondimen- 

sional perturbation terms of pressure, peripheral velocity, and axial flow rate. 

The cross-sectional area in the i^h gland is represented as 

ii,= i(#i + <Ti) (i5) 

where h± and 6± are the height of the gland and the radial labyrinth clear- 
ance . 


By denoting the displacement of the center of the rotor (x,y) as 


x = Ti c&dti. t> 
y jf = Yz Mw&V 

the area of the i c ^ chamber section is obtained as 


(16) 
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Because the change of state in the gland is isentropic change, the following 
relation is obtained: 


Therefore 
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From the above equations, the following linear equations are obtained: 
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Equations (19) and (20) are rewritten in the matrix form as 

+ f [ ncos Qt cosf + r*s,nQt sin if ) 

+ g ( n si n Qt cosf - r*cos.Qt sin if ) ( 21 ) 


where T, V, and B are (2k - 1) x (2k - 1) matrix and u, f and g are 
(2k - 1) row vectors and u is represented as 

llO = L , *1., $2 ,1* , — - , %*-l / 1i-l r 1 J 


ANALYSIS BY FINITE DIFFERENCE METHOD 
Multiplying equation (21) by T~1 yields 

I(U + Biu’ + (QlU= S ( Qt ok if + r 2 sinQt sin? ) 

+ r ( r, slnQtcss? - r 2 c»$Qt sin? ) < 22 ) 

where 

B = T"' V , G = T"’B 
S = T"JF , r =T~‘$ 

| = unit matrix 


By dividing the circumferential space of the rotor into n elements as shown 
in figure 3, the following finite difference equation is obtained about the j c 
element : 
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When the boundary conditions are set as u . (<(, j , t) = u, (<(>, + 2ir, t) , equation 
(23) is reduced for the overall system to J J 

IX +AX = dt cm S2 1> + lb Mm, £2 1 ■+ <c - dc^fS2.t> (24) 

where x, a, b, c, d are n(2k - 1) row vectors and A is n(2k-l) x n(2k-l) 
matrix 

The solution of equation (24) is obtained in matrix form as 

X-e**e + J' i e (t - s) UcA,J2sds+J o *e < *- s) \^J2sds 
+ J e ct S)A CMn&sds - f e J c#JJ2sds «5) 

y 0 J 0 

If the rotor is rotating at the steady state, the perturbation terms are equal 
to zero; so the initial condition for analysis becomes 

(£ = 0 (26) 
By using this condition, equation (25) becomes 

X =E r, gmSZi t Fr 2 + G (~r,£> ) 

t IH Y~zS2 (27) 

where 
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By using the nondimens ional variable of pressure, fluid forces acting on the 
rotor are given as 
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where E-jj is nondimensional variable of the i t ^ 1 gland and the j ^ element. 

By applying the relation 

r,c#fJ 2i> = x , r 2 x^vi2t = # (29 

= x , nJ2 cm S2t - 

to equation (28) , the flow- induced force due to the labyrinth seal is obtained 
as 
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These coefficients are the spring constants and damping coefficients for the 

gas flow through the labyrinth seal when the rotor moves parallel to the cas- 
mg axis. 


NUMERICAL EXAMPLE 

The labyrinth seal having three teeth is used as a numerical model, and 
the seal is divided into 24 elements for the finite difference method. Details 

! f u he , kby ! 1 ? th Sea1 ’ the rotor > and the fluid are shown in tables 1 to 3. In 
table 1, models A to C investigate the effect of seal clearance, models D to F 
investigate the effect of precession, and models G and H investigate the effect 
of divergence and convergence seals. 

For the calculation, the following values are used: 


Flow coefficient =0.7 


A = 64/Re for Re < 1200 
A = 0. 3164Re °* 25 f or Re > ^200 
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where Re is Reynolds number. 


Table 4 shows the pressure and flow rate in the steady state. Figures 4 
to 6 show the spring constants and damping coefficients for each model. Be 
cause the rotor axis is coincident with the casing axis for the calculation, 
the following relations are obtained: 

Kxx = Kyy, Kxy = -Kxy 

Cxx = Cyy, Cxy = -Cyx 

Figure 4 shows the spring constants and damping constants for models A, B, and 
C in order to investigate the effect of clearance. From figure (b) , pnncipa 
diagonal terms of the damping coefficients are constant for the variation of 
rotating speed; on the other hand, the cross terms of the damping coefficient 
vary with rotating speed. Decreasing the clearance makes the coefficients 
large. Figure 5 shows the effect of precession, from which it is known that 
the diagonal terms of both the spring and damping coefficients are constant for 
the variation of rotating speed but that the cross terms of both the spring and 
damping coefficients are strongly dependent on rotating speed. Figure 6 shows 
the different effects of convergence and divergence seals. From these figures 
it becomes clear that the cross terms of the spring coefficients and the damp 
ing coefficients are quantitatively different. 

ENERGETIC APPROACH TO STABILITY 


In this section, the energy, eigenvalue, and phase difference between two 
modes are derived for a two-degree-of-f reedom system, and the result is applie 
to the stability analysis of a rotor having a labyrinth seal. The general equa 
tion of motion is represented in matrix form as 


!Mic+Bi + IKx = ® 


(36) 


where 


1 = 


ri,i 

IM = 

"m u 

0 ' 

B = 

' b„ 

b* 

K = 

k H 

\ 

k 12 

IxJ 

/ 

.0 

m 2 * , 

I 

.b 2 , 

b 2 2. 

} 

s k 2 i 

k 22, 


As equation (36) has two eigenvalues, we consider the energy for each mode by 
assuming the periodical vibrational mode, that is, by neglecting t e rea P ar 
of the eigenvalue. Then the vibrational energy for one cycle for each mode 

written as 


E = f i T M idt + $ i T B id t + 5 


(37) 
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The first term of expression (37) is kinetic energy, the second is dissipative 
energy, and the third is potential energy. The energies of each term for one 
cycle (one period = t) are obtained as 

\ T i T M idt = o os) 

J 0 

K t B Kdt = ^ (b„x, 2 + b J2 i; +2bX,i 2 )dt = E, + e 2 (39) 

where 
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ft . r r 

\ * i^Xdt = \ Ak( X,x 2 - X,x 2 )dt = Ea ( 40 ) 
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where 

Ak — C ^2 ~~ K2 1 ) 

Therefore the total energy of the system for one cycle becomes 

E = E, + E 2 + E 3 , ( 4 D 

where 

E,= ^(^uf + b^uj), 

E 2 = ^*2 b u, u 2 c©s ( <£, - 4 > 2 ) t (42) 

E 3 = -2Ak a, u 2 s i n c 4 >, - <#> 2 ) , 

u i is the eigenvector and (<p^ - <p 2 ) is the phase angle between the first and 
second modes. By using this result, the stability of the system may be stated 
from the energy point of view as follows: 
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(1) If E > 0, the energy of the system is absorbed, and thus the system 
is stable. 

(2) If E < 0, the energy of the system is dissipative, and thus the sys- 
tem is unstable. 

In expression (42), E-^ is the energy obtained by diagonal elements of the 
damping coefficient. It is always positive; so if the damping coefficient is 
positive, this term always makes the system stable. Also E 2 is the energy 
obtained by cross elements of the damping coefficient, and its sign is depen- 
dent on the phase angle between <}>^ and 4*2 • Finally E 3 is the energy ob 
tained by cross elements of (k ^2 “ ^ 21 ^ anc * t ^ ie P^ ase angle between $1 and 42* 

From the above discussion, cross elements of the stiffness and the diagonal 
elements of the damping make the system unstable, and the cross elements of the 
damping do not affect the stability in this case. 

EXPERIMENTS 

The experimental apparatus shown in figure 7 was used to observe the flow 
pattern in the labyrinth seal and to investigate the dynamic behavior of the 
labyrinth seal. The rotor is driven by a variable-speed motor system, and its 
bearing (with eccentricity) is also driven by another variable-speed motor in 
order to obtain an arbitrary whirling speed and spinning speed. A two-stage 
labyrinth seal (straight type) is set up at the rotor, whose depth of gland, 
pitch, and mean clearance are 18.2, 30.0, and 1.8 millimeters, respectively. 
Rotating speed of the rotor and shaft for whirling drive are 84 337 and 
93 460 rpm, respectively. The casing is made from polymethyl-meta-acrylate in 
order to show the flow state, and water is used for the working fluid. 

Pressure in the shroud is measured by the semiconductor pressure gage, and 
its signal is analyzed by a real-time analyzer. Figure 8 shows the flow pat 
tern in the shroud, where a continuous vortex in the circumferential direction 
occurs in the fluid flow. The form is like a sinusoidal wave which is rotating 
in the same direction as the rotor. Figures 9 and 10 show the dynamic pressure 
(perturbation term) of the gland and the phase angle between deflection and the 
pressure for forward and backward precession, respectively. From these figures 
it is known that the dynamic pressure of the both cases is increased in propor 
tion to the increase of rotating speed but that the tendency of the phase angle 
of both cases is the reverse. 

MOTION PICTURE SUPPLEMENT 

The 8 -mm film was taken in order to observe the flow pattern in the gland 
of the labyrinth seal, where a continuous vortex in circumferential direction 
occurs. The vortex form is like a sinusoidal wave which is rotating in the 
same direction as the rotor. Figure 11 shows the two different flow patterns. 
The left side shows a conventional mathematical model of the flow pattern in the 
gland. This model is usually used to derive the fundamental equation. The real 
flow pattern is not similar to the model. In the gland the flow is composed of 
vortex and expansion flow as shown by the right side of figure 11 . 
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c Th j S the ? ext sut> j ect we should approach is the mathematical derivation of 
the fundamental equation considering the vortex in the gland. 

CONCLUSIONS 


The force induced by the labyrinth seal is solved by using the finite dif- 
ference method, expressed in proportional terms to the velocity and displace- 
ment of the rotor, and reduced to the expression of oil-film force of journal 
earings. Then the effect of the coefficients of the induced force on the 
stability of the rotor system is discussed from the energy point of view. The 
orce induced by the labyrinth seal always makes the rotor system unstable, and 
the tendency is remarkable when leakage of the seal is small. The flow pattern 
m the labyrinth seal was investigated experimentally, and it is known that a 
continuous vortex in the circumferential direction occurs in the fluid flow. 

The form is like a sinusoidal wave which is rotating in the same direction as 
the rotor. 


EDITORIAL SUPPLEMENT 

Reproduction of the film frame (fig. 8) was not very successful, so we 
introduce figures 12 and 13 in an attempt to demonstrate the film supplement. 
Figure 12 is a sketch of what appears in the film to be a spiral vortex. Fig- 
ure 13 represents a possible sequence of motions of the vortex center resulting 
from the periodic behavior of the flow interface. One must also be aware of 
the possible disturbance caused by the air bubbles in the flow field. 


These instabilities appear to be linked with those noted in unpublished 
work by Robert C. Hendricks and T. Trent Stetz of the NASA Lewis Research 
Center, where a flow visualization study was carried out on a water table to 
determine some characteristics of flows through sequential Borda-type inlets 
(no rotation and no centerbody). In figure 14(a) the four lucite Borda models 
were placed in such a way that they touched each other to form a continuous 
channel. The inlet water depth was similar to the passage width. The injected 
dye revealed that the flow through this configuration continued uninterrupted 
after passing the vena-contracta. The models were then placed with spacings of 
1/3 of the channel passage width (fig. 14(b)). The flow still continued in an 
uninterrupted manner after the vena-contracta. The models were then placed 

with spacings of 3/2 of the channel passage (fig. 14(c)). At this separation 
distance part of the flow entered the cavities and slight oscillations were ob- 
served. At a separation distance of 2-1/4 channel passage widths a very strong 
oscillation was observed (fig. 14(d)): The exhaust of one passage would "fan" 

the flow across the inlet of the subsequent Borda passage. These oscillations 
weakened when the separation distance was increased to 4 channel passage widths 
(fig. 14(e)). A separation of 6 channel widths showed minor oscillations (fig. 
14(f)). At distances beyond 16 channel widths the flow through each Borda pas- 
sage appeared nearly independent of the preceding flow (fig. 14(g)). 
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Figure 2. - Calculation procedure for steady state. 
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Figure 5. - Coefficients of flow- induced force for models D to F. 
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(a) Spring coefficients. 


(b) Damping coefficients. 


Figure 6. - Coefficients of flow-induced force for models G and H. 
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Figure 7. - Testing facility for labyrinth. 








Figure 8. - Flow pattern in gland. 
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Figure 9. - Experimental results for forward precession. 












(b) D/3 Separation - jetting. 





(e) 4D Separation, oscillation. 



(f) 6D Separation, minor oscillations. 



(g) 16D Separation, nearly independent of reservoir. 


Figure 14. - Flow visualization of four sequential Borda-type inlets. 
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